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Two Non-Hermitian fermion models are proposed and analyzed by using Foldy- 
Wouthuysen transformations. One model has Lorentz symmetry breaking and the other 
has a non-Hermitian mass term. It is shown that each model has real energies in a given 
region of parameter space, where they have a locally conserved current. 
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I. INTRODUCTION 

This work is concerned with fermion theories that are described by non-Hermitian Hamiltonian 
densities. (For earlier studies, see Refs. [ll-[^.) We are particularly interested in fermion models 
that violate Lorentz invariance. 

Our general procedure is to begin with a non-Hermitian Hamiltonian density of the form H = 
for which the Hamiltonian operator PL has eigenvalues u, which are assumed to be real in a 
region of parameter space. The Schrodinger form of the equation of motion is 

idotp = T-Lip. 

We then introduce a mapping U to implement the Foldy-Wouthuysen (FW) transformation on 
the free fermion field ^|J and PL, such that 

X = Uif and UPLU~^ = 

The equation of motion can then be rewritten as 

idoX = W7°x, 

where the FW Hamiltonian density is PIfw = = ^XX- 

The motivation for this construction is to map the non-Hermitian Hamiltonian density H onto 
the FW Hamiltonian density Hfw^ which is Hermitian (for real ui), with the same eigenvalues uj 
as the original Hamiltonian operator PL. Unlike the Dirac case, the operator U in this paper is not 
unitary, and thus Hfw 7^ H. Indeed, 

H = = iipw- ( 1 ) 

In Ref. it is shown that the FW transformation for an extension of the Dirac equation that 
contains all the Lorentz-symmetry and CRT-violating terms is consistent with the Standard Model 
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Extension (SME) [^. However, Ref. deals with Hermitian Lagrangians and does not overlap 
with the present stndy. We also note that a non-Hermitian Lagrangian is studied in Ref. [^, for 
which the corresponding Hamiltonian is Hermitian. This feature is possible because of the presence 
of mixed derivatives dodk in the Lagrangian. 


A. Derivation of the equation of motion 


Before proceeding, it is important to explain the derivation of the equation of motion. We will 
show that the equation of motion may be obtained by the usual (formal) procedure of performing a 
variational derivative of the action with respect to ^ and not varying ip. (Of course, this procedure 
is questionable for a non-Hermitian Lagrangian because varying the action with respect to ip does 
not reproduce the same equation of motion after taking the Hermitian conjugate.) 

The four fermion components may be written as ipa = 4>a + iXa, where cpa and Xa are real. The 
action then has the form 

S = J Ip {i0 - m - r) Ip = J {(pb - iXb) {4>c + iXc), 


where T satisfies Tf = T and { 7 °,r} = 0. 

The equations of motion can then be correctly obtained by demanding that the variational 
derivatives of S with respect to (pa and Xa independently vanish. We obtain 


5S 

S(pa 


- 7 °r)^^ {(p + ix)c - - ^ 7 ° - 7 °r)^^ {(p - ix)b 

i7V + i ( 7 °^) ^ - m7° + m7°^ - 7°r + (7°r)^j 




— irri'y^ — — i7'^r — i(7°r) 


J ac 
0r\T 


Xc- 


Then, from the properties of gamma matrices in the Dirac representation, we get 


6S 

S(pa 


2 

+2 


*7°(7°5o + 7^^i + 7^^3) - 7°r 


1 + e 


—7^7^52 — imx^ — i7'^r 


0 „-„ 0 t 


1 - e 


Xc, 


where e = -1-1 if T is real and e = — 1 if T is imaginary. Similar steps lead to 


bS 

SXa 


-2i 

+2 


^7^7^52 — — 7^r 


0 „, 0 t 


1 - e' 


i7°(7°5o + X^di + 7^53) - 7°r 


1 + e 


Xc- 


It is now easy to see that 


1 bS i bS 

2b(pa~^ 2 bXa 


-m-r){(p + ix))^, 


which is independent of e. The latter equation is equivalent to 7^^- = ( 1 $ — m — T)^p, which 
agrees with the variation bS/bip performed with ^p and ip considered as independent fields. Hence, 
we may derive the equations of motion by following the formal conventional procedure. 
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II. MODEL OF LORENTZ-SYMMETRY VIOLATION 


In this section we study a specific form of Lorentz-symmetry-violating kinematics for (3 + 
l)-dimensional fermions, which have real energies if the parameters of the model satisfy certain 
inequalities. Using the metric = diag(l, —1, —1, —1), we begin with the Lagrangian 


C = ijj (i0 — — m) 


( 2 ) 


where 6^ represents the (real) vacuum expectation value of a vector and this signals the breaking 
of Lorentz symmetry. The term (without the factor of i) appears naturally in the SME. 
This term can be obtained from the QED one-loop effective action in curved spacetime j^, where 
b,, oc Oi is the fine-structure constant, m is the electron mass, and R is the Ricci scalar for 

the curved background. 

The Hamiltonian density corresponding to the Lagrangian ([2|) is 


H = —-ip — C = pj [ij ■ V + ipi + m) 'tp, 
dip ^ ' 

where ip = d^ip. Note that H is not Hermitian because it contains the anti-Hermitian term 
iipjpip = — (iippiipy. The dispersion relation corresponding to the Lagrangian ([2]) is 

{u! — = m? + {p — ib)"^. 

However, if we choose 5o = 0 and restrict our attention to motion in the plane perpendicular to b, 
the energies uj satisfy 


uj‘^ = m? + p^ — b^ for 6o = 0 = p • 6, 


(3) 


and they are real for any momentum p as long as b^ < vn?. The existence of imaginary energies 
for low momentum (when 6^ > mp] implies the possibility of runaway modes, which is a problem 
known in SME studies (see Ref. 10| for a recent discussion). (Runaway modes in RT-symmetric 
theories are also considered in Ref. |lll|.) 


A. Current conservation 

To verify current conservation we multiply the equation of motion on the left by ip and then 
multiply the Hermitian conjugate of the equation of motion on the right by 'y^ip gives 

iip{0 — b ■ ^)ip = mipip and — iip{ 0 +b ■ j)ip = mipip. 

Next, we subtract these two equations and obtain 

dfj_{ip^^ip) + 2ip{b-^)ip = 0. 

As usual, we define the probability density as p = ip'y^ip = ip^ip and the current density as 
j = 'ipjip. The previous equation now reads 

p-V-j = -2b-j. (4) 

Because we restrict our attention to motion in the plane perpendicular to b, the current j is 
perpendicular to b and thus b ■ j = 0. Therefore, the continuity equation dU has the usual form 
d^{ip'y^ip) = 0 and the current ip^^^ip is conserved. 
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B. Charge-conjugation, parity, and time-reversal symmetries 


In the framework of the SME, the Lorentz-symmetry-violating term ^ 6-7 t/j (without the factor 
of i), is CVT-odd 1^. In our case, with b —>■ ib, we can check each symmetry independently: 

• Charge-conjugation C: 

b ■ 7 b ■ 7 "^ C)^l: = —5 • 7 V’j 


, such that Ip ib ■ ^ Ip is C-odd. 


where C satisfies C ^^fiC = ^ = — 7 ^ 

• Parity V: 

'’~^{i'ip b ■ 7 'ip)V = ^' 07 ^ b • 7 'y^'ip = —ip) 6 • 7 ' 0 , 


V- 


such that Ip ib ■ j ip is "P-odd. 

• Time-reversal T : 

7 ^{ip) b ■ 7 p))T = —ip)'y^C~^ b ■ ^ C'j^p) = ip)'^^ b ■ 7 ^ 7 ®^ = i'lp b ■ p p), 
such that p) ib ■ p p) is 7 ~-even. 

The non-Hermitian term ip) b ■ p p) is therefore even under CVP but odd under VT. 


C. Foldy-Wouthuysen transformation 


The Schrodinger form of the equation of motion obtained from the Lagrangian ([2]) is 


idop) = Tip), 

where the Hamiltonian operator is P = 7 ^ {p — ib)-p + m 
consists of writing the equation of motion in the form 

'idox{t,p) = u)x^x{t,p), 


The FW field transformation 


(5) 


where u) is the energy obtained from the dispersion relation (l3|), and y = Up), where U is to be 
determined. This form of the equation of motion explicitly shows the evolution of the positive and 
negative energy modes, and leads to the Hamiltonian density 

= X^W7°x = U)XX, 


which is Hermitian as long as u) is real. We will use the FW transformation to map the non- 
Hermitian Hamiltonian density H, with eigenvalues u), to the Hermitian Hamiltonian HpW) with 
the same eigenvalues. 

Since the Hamiltonian density TL has a similar structure as in the Dirac case, by analogy with 
the latter, we seek U in the form 


U = exp 


{p- ib) ■ p 

y/p2 _ 52 


U-ib) ■ p 

= cos 9 + ^ sin 9, 

_ 52 


( 6 ) 


where is imaginary in the infrared regime < 6 ^. Unlike the Dirac case, U is not 

unitary: U\9) 7 ^ U{—9) = U~^{9). In the ansatz ([ 6 j), the FW angle 9 must be determined in order 
to obtain the equation of motion ([5]). It is then easy to see that 


unu-^ 


7 °(p -Tib) -p 


cos( 20 ) — 


m 


Y^p2 — 52 


sin( 20 ) 


+7 


mcos(29) + \/p^ — b^ sin(20) 
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Canceling the term proportional to 7 requires that 

tan(26*) = y/p"^ — j m. 

For 0 < Re (20) < 7 r /2 this implies that 

cos( 20 ) = m/w and sin( 20 ) = sj p^ — juj. 

This result reduces to the usual definition of the FW angle tan(20) = p/m in the Dirac case, and 
in the present case it is then straightforward to check that one obtains the expected form 


UT-LU ^= 7 ° mcos(20) + ^/p^ — sin(20) 


= W7°. 


It is interesting that even when iP' < m? the infrared regime p^ < iP is characterized by a purely 
imaginary FW angle although the energies oj are real. The transition between the real and purely 
imaginary FW angle corresponds to the limit where 0 —0. 


III. PARITY-VIOLATING MASS TERM 

Let us examine the FW transformation for a model that was looked at in the context of PT- 
symmetric theories [H, and whose features are similar to those of the model examined in Sec. m 
although there was no Lorentz-symmetry violation. The Lagrangian considered in Ref. [l| is 

C = 'ip {i$ — m — ip, (7) 

which contains the anti-Hermitian mass term pip'j^'ip = —p{Pp'^^'ip)P One can check that the 
corresponding Hamiltonian density is 

H = •ip • V + m + P 7 ^^ ip, 

and is RT-even, but R-odd and T-odd. While the Hamiltonian is not Hermitian, the dispersion 
relation 


2 2,2 2 /o^ 

(jj = m + p — p (8) 

shows that the energies are real for any momentum as long as mP > pP. The dispersion relation 
(jHjl is the same as that for ([3]), when we substitute b'^ ^ pP. As a consequence, we expect the FW 
transformation to be similar to that described in the previous section. 

The Schrodinger form of the equation of motion obtained from the Lagrangian ([7]) is 


idoip = PLip, (9) 

where the Hamiltonian operator is 77 = 7 ® [p • 7 + m + P 7 ^]. By analogy with the model ([2]), we 
look for a mapping of the form 


U = exp 


p-p + PY 

^/jP^PpP 


= cos 0 + 


p • 7 + /i7 


sin 0 . 


y/jpPPpP 

The parameter 0 should be determined in order to write the equation of motion in the form 

'idox{t,p) = ujj°x{t,p), 


( 10 ) 


( 11 ) 
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where x^Uxjj and iv satisfies the dispersion relation ([ 8 ]). Analysis similar to that described in the 
previous section shows that one must choose 

tan( 20 ) = \/p^ — IJ? jm. ( 12 ) 

Hence, the operator (IlOp maps the non-Hermitian Hamiltonian density H to the FW Hamiltonian 
density Hpw = loxx^ which is Hermitian as long as > fj?, and has the same eigenvalues ui 
as H. (The present FW mapping is different from that described in [l|, where U contained the 
non-Hermitian term only, and therefore led to a different Hermitian Hamiltonian density.) 


A. Alternative description 


Even in the situation m? > for which the energies are real, one can distinguish the UV 
regime > p?, where the angle 9 defined by (jl2p is real, from the IR regime < p?, where 9 is 
purely imaginary. The identity (fToP is based on the property that (p • 7 + P 7 ^)^ = — p^. One 

may consider instead 


U' 


exp 


(q, P^jJrpr^ 
\ _ p2 j 


cosh 9' + 


\/p^ -P^ 


for which the parameter 9' must be chosen so that 

tanh( 20 ') = ■\/ p^ — p^/m 


(13) 


in order to obtain the FW form m of the equation of motion. One can easily see that the 
mappings m and (jl3p are equivalent if one identifies 9 = i9', which corresponds to changing the 
IR regime to the UV regime < p^ in the FW mapping. 

Finally, we note that 


lim U = lim U' 


^ p-7 + P7^ 

2m 


so that the apparent singularity when p^ —>■ p^ is actually a smooth limit, as one would expect 
from physical grounds, and the FW transformation is not singular at p^ = p^. 


B. Conserved current 


The eigenstate solutions of the equation of motion Q satisfy 

'^{t^P) = exp(—i?^t)V'( 0 ,p). 


Thus, 


'ip^{t,p) = V'^( 0 ,p) exp(f'R^t), 

and the corresponding naive probability density does not respect unitarity because 
'il)\t,p)'il){t,p) ='!/’^( 0 ,p) exp(fHl’f) exp(—z?{f)V'( 0 ,p) 7 ^ V''^( 0 ,p)V'( 0 ,p). 

For consistency, we therefore need to redefine the probability density, which we write in the form 

p{t) = V^^(t)AV^(t), 
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where A is to be determined so that p(t) = p(0) for eigenstates. This leads to 

Texp(—= exp(—i9{h)A. 

The latter identity implies that ATi = T-0A, which is solved by expanding A in the basis ( 1 , 7 , 7 ®). 
We find that ^ = 1 + ^ 7 ®, and the probability density, which respects unitarity, is 

V m / 

From this result we find the conserved current as follows. For a general solution to the equation 
of motion ([9]) we have 

p = P (l — 7 ®') -(/; + 7 + — 7 ^) V’, 

and according to the same equation of motion 

= V • 7 + — /r 7 ®) and -i/) = 7^7 • V'ljj — 

As a consequence, we see that /? = V • [■0 7(1 + ^ 7 ®) V'] • Thus, the conserved current is 

f = (1 + — 7 ^) i’- 

Finally, we decompose 'll: into right- and left-handed components, xf: = 'i/’R + V'L, by using the 
projection operators ^ (l ± 7 ®): 

V’R = Hi + H) V’, V’L = Hi “ H) V'- 

Note that V’R V'L are eigenstates of 7 ® ( 7 ®V'r = V’R; 7 ^V’L = —'4’l) and that they are 
orthogonal (V’rHl = V’lHr = 0 ). We then observe that the probability density becomes 

P = (1 + —) IV’Rp + (1 “ —) iV'Lp, 

which is always positive in the relevant regime where < m?. Thus, p = m and p = —m 
are interesting special cases in which the contribution to the density is entirely from right- or 
left-handed degrees of freedom. 

To conclude, we have discussed two potentially relevant applications of the non-Hermitian mass 
term /ry®. First, because right and left-handed helicities do not contribute in the same way to the 
probability density, this model could shed new light on neutrino physics. A dynamical mechanism 
to generate this non-Hermitian mass term is one possible extension of this work. It would then be 
interesting to gauge this model and study the effect of the non-Hermitian mass term on the chiral 
anomaly, which might be cancelled independently of the number of lepton generations. 
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